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Abstract 

This study proposes an approach based on a perturbation technique to con¬ 
struct global solutions to dynamic stochastic general equilibrium models 
(DSGE). The main idea is to expand a solution in a series of powers of a 
small parameter scaling the uncertainty in the economy around a solution 
to the deterministic model, i.e. the model where the volatility of the shocks 
vanishes. If a deterministic path is global in state variables, then so are the 
constructed solutions to the stochastic model, whereas these solutions are 
local in the scaling parameter. Under the assumption that a deterministic 
path is already known the higher order terms in the expansion are obtained 
recursively by solving linear rational expectations models with time-varying 
parameters. The present work also proposes a method rested on backward 
recursion for solving general systems of linear rational expectations models 
with time-varying parameters and determines the conditions under which the 
solutions of the method exist. 



1 Introduction 


Perturbation methods applied in macroeconomics are used to expand the ex¬ 
act solution around a deterministic steady state in powers of state variables 
and a small parameter scaling the uncertainty in the economy. The solu¬ 
tions based on the Taylor series expansion are intrinsically local, i.e. they 
are accurate in some neighborhood (presumably small) of the determinis¬ 
tic steady state. Out of the neighborhood the solutions may behave odd, 
for example, can imply explosive dynamics (Kim et ah , 2008). The other 
problem with the perturbation method is that we do not know a priori for 
non-trivial models how small the neighborhood must be to achieve a given 
level of accuracy. 

The recent crisis has renewed an interest in methods that provide global 
solutions to DSGE models, i.e. the solutions some points of which are far 
away from the steady state. This may occur after a big shock hitting the 
economy, or if the initial conditions are far away from the steady state, the 
examples of this situation are the economies in transition and developing 
economies. 

This study presents an approach based on a perturbation technique to 
construct global solutions to DSGE models. The proposed solutions are 
represented as a series in powers of a small parameter a scaling the covariance 
matrix of the shocks. The zero order approximation corresponds to the 
solution to the deterministic model, because all shocks vanish as cr = 0. 
Global solutions to deterministic models can be obtained reasonably fast by 
effective numerical methods^ even for large size models (Hollinger , 2008). 
For this reason the next stages of the method are implemented assuming 
that the solution to the deterministic model under given initial conditions is 
known. 

Higher-order systems depend only on quantities of lower orders, hence 
they can be solved recursively. The homogeneous part of these systems is the 
same for all orders and depends on the deterministic solution. Gonsequently, 
each system can be represented as a rational expectation model with time- 
varying parameters. In the case of rational expectations models with constant 
parameters the stable block of equations can be isolated and solved forward. 

^The algorithms incorporated in the widely-used software such as Dynare (and less 
available Troll) find a stacked-time solution and are based on Newton’s method combined 
with sparse-matrix techniques (Adjemian, Bastani, Juillard, Karame, Mihoubi, Perendial, 
Pfeifer, Ratto, and Villemot , 2011). 
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This is not possible for models with time-varying parameters. 

The other contribution of the present work is a method proposed for solv¬ 
ing general systems of linear rational expectations models with time-varying 
parameters and determines the conditions under which the solutions of the 
method exist. The method starts with finding a hnite-horizon solution by 
using backward recursion. Next we prove that as the horizon tends to inhn- 
ity the finite-horizon solutions approach to a limit solution that is bounded 
for all positive time. The proposed method for solving linear rational expec¬ 
tations models with time-varying parameters may be valuable in itself, for 
example, for solving models with anticipated structural and policy changes. 

Notice that whenever the deterministic solution is global in state variables 
so is the approximate solution to the stochastic problem. At the same time, 
if the parameter a is small enough, then the solution obtained is close to the 
deterministic one. For this reason, we shall call this approach semi-global. 

To illustrating how the method works we apply it to the asset pricing 
model of Burnside (1998). The simplicity of the model allows for obtaining 
the approximations in an analytical form. We compare the policy functions 
of the second order solution of the semi-global method with the local Taylor 
series expansion of orders two and six (Schmitt-Grohe, and Uribe , 2004). 
The results show that the semi-global solution is more accurate, in some 
sense, than even the sixth order of the local Taylor expansion. 

This paper contributes to a growing literature on using the perturba¬ 
tion technique for solving DSGE models. The perturbation methodology in 
economics has been advanced by Judd and co-authors as in Judd (1998); 
Gaspar, and Judd (1997); Judd, and Gnu (1997). Jin, and Judd (2002) 
give a theoretical basis for using perturbation methods in DSGE modeling; 
namely, applying the implicit function theorem, they prove that the per¬ 
turbed rational expectations solution continuously depends on a parameter 
and therefore tends to the deterministic solution as the parameter tends to 
zero. 

Almost all of the literature is concerned with the approximations around 
the steady state as in Gollard, and Juillard (2001); Schmitt-Grohe, and 
Uribe (2004); Kim et ah (2008); Gomme, and Klein (2011). Lombardo 
(2010) and Lombardo, and Uhlig (2014) make use of series expansion in 
powers of a to provide a theoretical foundation for pruning methods (Kim et 
ah , 2008), which is aimed to avoid the explosive behavior of a solution. An- 
dreasen, Eernandez-Villaverde, and Rubio-Ramrez (2013) develop the same 
approach for higher-order approximations. Lombardo and Uhlig’s approach 
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can be treated as a special case of the method proposed in the current study, 
namely a deterministic solution around which the expansion is used is only 
the steady state. Both approaches based on the perturbation methodology 
used in applied mathematics ( Nayfeh (1973) and Holmes (2013)). The 
essence of the methodology is to expand a solution in a series of powers of 
a small parameter, and thus obtain a set of problems that can be solved 
recursively. It is supposed that each of this problems is easier to solve than 
the original one. Actually in applied mathematics literature (Nayfeh (1973) 
and Holmes (2013)), the zeroth-order approximation is typically a function 
of time t rather than a steady state as in Lombardo (2010), Lombardo, 
and Uhlig (2014) and Andreasen, Fernandez-Villaverde, and Rubio-Ramrez 
(2013). Judd (1998, Chapter 13) outlines how to apply perturbations around 
the known entire solution, which is not necessarily the steady state. He con¬ 
siders a simple continuous-time stochastic growth models in the dynamic 
programming framework. This paper develops an approach to construct ap¬ 
proximate solutions to discrete-time DSGE models in general form by using 
the perturbation method around a global deterministic path. 

Despite the fact that the pruning procedure avoids the explosive behavior 
of a solution, it remains local, and as such may have some undesirable prop¬ 
erties. For example, the pruning procedure might provide a first few impulse 
responses with wrong signs under a sufficiently large shock. This case seems 
even worse than the explosive dynamics since the impulse responses for a first 
few periods are most interesting and relevant for theoretical implications of a 
model as well as a policy analysis; therefore, their incorrect signs could mis¬ 
lead a researcher or a policymaker. In this situation the pruning procedure 
just conceals the real problem. As we will show in the example, the problem 
with a wrong sign of impulse responses can occur even in a situation where 
the pruning is not needed. 

The rest of the paper is organized as follows. The next section presents the 
model set-up. Section 3 provides a detailed exposition of series expansions for 
DSGE models. In Section 4 we transform the model into a convenient form 
to deal with. Section 5 presents the method for solving rational expectations 
models for time-varying parameters. The proposed method is applied to an 
asset pricing model in Section 6, where it is also compared with the local 
Taylor series expansions. Gonclusions are presented in Section 7. 
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2 The Model 

DSGE models usually have the form 


Etf{yt+i, yt, xt+i, xt, zt+i,zt) = 0, 

Zt+l = Ezt + cr£t+i, 


( 1 ) 

( 2 ) 


where Et denotes the conditional expectations operator, xt is an Ua, x 1 vector 
containing the t-period endogenous state variables; yt is an x 1 vector 
containing the t-period endogenous variables that are not state variables; Zt 
is an X 1 vector containing the t-period exogenous state variables; St is the 
vector with the corresponding innovations; and the x covariance matrix 
fl; / maps x x x x x into x and is assumed 
to be sufficiently smooth. The scalar cr (cr > 0) is a scaling parameter for 
the disturbance terms et- We assume that all mixed moments of St are hnite. 
All eigenvalues of the matrix A have modulus less than one. The problem is 
to hnd a stable solution {xt, yt) to (1) for a given initial condition (xq, Zq). A 
process is stable if its unconditional expectations are bounded (Klein , 2000). 

3 Series Expansion 

3.1 The General Case 

In this section we shall follow the perturbation methodology used in applied 
mathematics (see, for example, Nayfeh (1973) and Holmes (2013)) to derive 
an approximate solution to the model (1)~(2). For small a, we assume that 
the solution^ has a particular form of expansions 



(3) 



(4) 


The exogenous process Zt can also easily be represented in the form of ex¬ 
pansion in a 



(5) 


^As is conventional in applied mathematics literature (see for example Holmes (2013)), 
each term of the expansion is not divided by the factorial term at this stage. 
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Indeed, plugging (5) into (2) gives 

zt+i = = ^( 4 °^ + 


Collecting the terms of like powers of a and equating them to zero, we get 

44i = (6) 

4+1 = ^4 ^ + ^i+i- (7) 

Since the expansion (5) must be valid for all a at the initial time t = 0, the 
initial conditions are 

4°^ = Zq and 4^^ = 0. (8) 

It is worth noting that in contrast to Lombardo, and Uhlig (2014), where 
= y and xf^ = x are the steady state, here xf''^ and y^^ are functions of 
t and, as will be shown below, are a deterministic solution of the problem. 
In this way the pruning procedure is always local around the steady state, 
whereas we focus on the case where and y^^ are global. 

Then substituting (3), (4) and (5) into (1), we have 

EtfU+i + a44 + (T^yi% + ■■■, + • • • , 

4+\ + + • • •, 4°^ + (^4^^ + + • • •, (9) 

4?i + (^z44, 4°^ + (^z4^) = 0. 

Expanding the left hand side of (9) for small cr, collecting the terms of 
like powers of a and setting their coefficients to zero, we obtain 

Coefficient of <7° 


,,(0) .^(0) .^(0) .(0) (0)x _ . 


( 10 ) 


The requirement that (4) and (5) must hold for all arbitrary small a implies 
that the initial conditions for (10) are 

4°^ = Zq and x® = Xq. (11) 

The terminal conditions for y^ and are the deterministic steady states 

2/^^ = y and = x. (12) 


The system of equations (6) and (10) is a deterministic model since it cor¬ 
responds to the model (1) and (2), where all shocks vanish (for this reason 
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we omit the expectations operator in (10)). The deterministic model (6) 
and (10) with the initial and terminal conditions (11) and (12), respectively, 
can be solved globally by a number of effective algorithms, for example the 
extended path method ( Fair, and Taylor (1983)) or a Newton-like method 
(for example, Juillard (1996)). As this study is primarily concerned with 
stochastic models, in what follows we suppose that the deterministic model 
is already solved and its solution is known. 

Coefficient of 

■ yf+i + /2,t • yf'’ + . 

+/5,t4i-i + = 0- 


The matrices 


fi,t = fi (z/m, yi 


(0) .^(0) ^(0) ^(0) ^(0) 


' ^t+ll ^ 




. i=l,. 


6 , 


are the Jacobian matrices of the mapping / with respect to the ith argu¬ 
ment (that is yt+i, yt, Xt+i, Xt, Zt+i, and Zt, respectively), at the point 

(^2/t+i) 1/1°^ 4+1, 4°^) • The requirement that (4) and (5) must hold 

for all arbitrary small a implies that the initial condition for (13) is 


Zq^'^ = 0 and = 0. 


(14) 


Coefficient of n > 1 

■ 4?i + /2,t ■ yt^ + f3,t ■ + U,t ■ + 4+U = 0. (15) 


The requirement that (4) must hold for all arbitrary small a implies that the 
initial condition for (15) is 

= 0. (16) 

A nice feature of the set of systems of equations (15) is that the linear 
homogeneous part fi^t is the same for all n > 0. The difference is only in the 
non-homogeneous terms fh^it are some mappings for which the set of 

arguments includes only quantities of order less than n 


,,( 0 ) ( 0 ) ( 0 ) ( 0 ) 
yt+iy yt ! 


7 ( 0 ) ^( 0 ) ^( 1 ) „(i) 

I yt+i y yt )-^t+i ! 


Particularly, for n = 1, 2 we have 

= (/wA + 
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and 


,{i) 

t+i 


Etvi% = Et{lfn,t (l/m) + \f22,t + \hz,t (a 

+ ^/44,i + lf55,t + lfe6,t + fi2,tyi+iyl^^ 

+fi3,tyi+i4+i + + fi5A+A+i + fw,tyi+i4^^ 

+ f23,tyi^^4% + f24,tyt^^x[^^ + f25,tyi^^4+l + f26,ty4^z4^ 

+ f34,txl^+ixi^'^ + f35,tx[^+l4+l + f36,tx[^+iz4'^ + f45 

+ f46,txi^^z4^ + /56,t4+l4^^}) 


( 17 ) 


respectively; where fij^t, * = 1, ■ ■ ■, 6, j = 1,..., 6, denotes the mixed partial 
Frechet derivative of ft of order two with respect to zth and jth argnments 
at the point 


,,{ 0 ) ( 0 ) ( 0 ) ( 0 ) ( 0 ) ( 0 ) 
yt+11 yt iXt+iiXt : ^t+1 ) 


(18) 


In other words, fij^t is a bilinear mapping (see, for example, Abraham, Mars- 
den, and Ratin (2001, p. 55)) depending on vector (18) (and hence on t).^ 
The expectations Etrjt+i d^ie bonnded if all conditional mixed moments of 
z44i are bonnded np to order n and the vectors (18) are bonnded for all 
t > 0. 

Eqnation (15) with the initial conditions (16) is a linear rational expecta¬ 
tions model with time-varying parameters and bonnded the non-homogeneons 
terms Ety^^i- To solve the problem (15)-(14) is eqnivalent to Ending a 
bonnded solntion {x4\y4'^) for t > 0 nnder the assnmption that the bonnded 
solntions to the problems of all orders less than n are already known. Know¬ 
ing how to solve these types of model and nsing the structnre of mappings 
Ety^^ii recnrsively solntions, y4\x4\ to (15) for every order n, 

starting with n = 1. 


3.2 An Example of the Series Expansion: An Asset 
Pricing Model 

In this section the method of expansion aronnd a deterministic path applies 
to a simple nonlinear asset pricing model proposed by Bnrnside (1998) and 
analyzed by Collard, and Jnillard (2001); Schmitt-Grohe, and Uribe (2004). 

^We do not make use tensor notation for brevity. 
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In this model the representative agent maximizes the lifetime utility function 


max 



Cl 

9 


subject to 


PtCt+i + Ct= PtCt + dtCt, 


where /? > 0 is a subjective discount factor, 6 < 1 and 6 ^ 0, Ct denotes 
consumption, pt is the price at date f of a unit of the asset, et represents 
units of a single asset held at the beginning of period f, and dt is dividends 
per asset in period t. The growth of rate of the dividends follows an AR(1) 
process‘d 

xt = (1- p) X + pxt-i + aet+i, (19) 

where Xt = ln{dt/dt-i), and et+i ~ NIID{0,1)- The hrst order condition 
and market clearing yields the equilibrium condition 


yt = f3Et [exp(6'a:t+i) (1 + yt+i )], (20) 

where yt = Pt/dt is the price-dividend ratio. This equation has an exact 
solution of the form (Burnside , 1998) 

OO 

yt = '^/d' exp [tti + bi{xt - S)], (21) 

i=l 

where 

1 / V 

at = 9x1 + - [ - - 

2\l-pJ 

and 

i-p 

It follows from (20) that the deterministic steady state of the economy is 

(3exp{9x) 

^ ~ 1- ^exp{9x)' 

^By abuse of our previous notation, we let Xt stand for the exogenous process as in 
Burnside (1998); Collard, and Juillard (2001); Schmitt-Grohe, and Uribe (2004). 


2p(l 


P^) I P^(l 


P 




P 


P^ 


( 22 ) 










We now express a solution to the system (19)-(20) as an expansion in pow¬ 
ers of the parameter a up to a second-order approximation and decompose 
the original problem into a set of auxiliary problems. Specihcally, assume 
that the solution can be represented in the form: 


(0) I (1) I 2 (2) 

yt = yt + m +(^ yl 

Xt = 


(23) 

(24) 


Substituting (24) into (19) and collecting the terms containing ct° and cr^, we 
obtain the representation (24) for xt 


X 


( 0 ) 

t+i 


= (1- p) X + pxl 


( 0 ) 


( 1 ) ( 1 ) , 
xl+i = pxy + Et+l. 


(25) 

(26) 


Since the expansion (24) must be valid for all a at the initial time t = 0, the 
initial conditions are 

Xq^^ = xo and = 0. (27) 

Substituting now (23) and (24) into (20) yields 


( 0 ) , ( 1 ) , 2 ( 2 ) , 

yt +(^ yi + 


= l3Et{e^p 


9 ( + ax 


( 1 ) 

t+i 


[1 + yf'’ + + ■■■]} 


,(i) 


2„(2) 


Expanding exponential for small a gives 


(0) I (1) I 2 (2) , 

yt +(^yt +(^ yl + ■■■ 


= (3Et exp{9x 


( 0 ) ^ 
t+i) 


1 + ^ [(^9x\:^-^j +■■■ 




[1 + 2/t+i + (Tyl^+i + cr'^yl+i + ■ ■ ■ ] 


Collecting the terms of like powers of a in the last equation, we have 

Coefficient of a^ 


yf^ = P exp(Px)+\)(l + ?/)+\), 


( 0 ) 


,(o) 


.^(0) _ _(o) 

Xt+l — PXt ■ 


(28) 

(29) 


Coefficient of a^ 

^ = exp(PxS+\)/?P (l + -F exp{ex[%)/3Etyl% (30) 
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( 1 ) ( 1 ) , 

xJVi = p^t + ^t+1- 


Coefficient of 

y? = |/5exp(0xJ+\)02 (l + Et (^xS+\) 

+6/3 exp{ex[%)Et (xS+\2/J+\) + /3 exp{ex[%)Et (yi%) . 


( 31 ) 


(32) 


The system (28) and (29) is a deterministic model. Its solution can easily 
be obtained by, for example, forward induction 


( 0 ) 

Vt = 


2 = 1 


exp 


6 


, p(i 

XI H-;- [Xt — X) 


P 


(33) 


Under the assumption that and x/^^ are already known for t ^ 0, 
Equations (30) and (31) constitute a linear rational expectations model with 
time varying deterministic coefficients exp(Px|°\)/3. The expectations of the 
term in (15) has the form exp{dx[^i)/3Et 6*xj)j.\(l + ?/|°\) . Equation 
(32) is also a linear forward-looking equation with time varying deterministic 
coefficients exp(6*x[°\)/d, and the term 


( 0 ) 


= /?exp(Px)+\)P 


(0) 


^6 (^1 + Et + Et 


depending only on solutions of orders less than two, i.e. yl+i, Pt+i- 
Therefore, both the system (30) and (31), and Equation (32) are linear 
forward-looking models with time varying coefficients. Under the condition 
that we know how to solve these types of model, they can be solved recur¬ 
sively starting with solving (30) and (31), then passing to (32). In Section 5 
we present a method for solving such types of model and prove the conver¬ 
gence of the solutions implied by the method to the exact solution. In the 
next section we transform equation (15) in a more convenient form to deal 
with. 


4 Transformation of the Model 

Define the deterministic steady state as vectors {y, x, 0) such that 

f{y,y,x,x,0,0) = 0. (34) 
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We can represent in (15) as fi^t = fi + fi,t, i = 1, • • •, 6, where 


fi = fi{y,y,x,x,0,0) 


are the Jacobian matrices of the mapping / at the steady state with respect 
to ith. argument, and 



ktivt+i, yi 


( 0 ) ..( 0 ) ^( 0 ) ^( 0 ) ^( 0 ) ^( 0 ) 


' ^t+ii ^ 


t+ii 


fi{y, V, x, x, 0,0). 


(35) 


Note also that fi^t —)• 0 as t —>■ oo, because a deterministic solution must tend 
to the deterministic steady state as t tends to inhnity. Consequently, fi^t can 
be thought of as a perturbation of /j. As Equations (15) have the same form 
for all n > 0, to shorten notation, further on we omit the superscript (n) 
when no confusion can arise. Therefore Equations (15) can be written in the 
vector form 


^tEt 


Xt+l 

yt+i 



+ Etr]t+i, 


(36) 


where = 


fs + /s.t, /i + /i,t 


fi + fi,t, f2 + f2,t 


We assume 


and At = 

that the matrices <ht are invertible for all f > 0. For instance, this assump¬ 
tion always holds in some neighborhood of the steady state if the Jacobian 
[/ 3 )/i]~^ the steady state is invertible.^ 

Pre-multiplying (36) by we get 


Et 


Xt+l 

yt+i 



+ Mt 


Xt 

yt 


+ ^t^Etrit+i, 


(37) 


where L = [/g, /i] ^ [U, /s] and 


Mt = 


-1 r 


/a + /s.t, fl + fl,t fi + fi,t, f2 + f2,t — [/s, fl] ^ [fi, f2] ■ 


Particularly, for n = 1 we have 


Et 

Ai) 

Xt+l 
. (1) 

= L 

Xt 

+ Mt 

^(1) 

Xt 

. (b 


yt+i 


[ yt \ 


[ yt \ 


+ ^(/5,f^ + /6,t)4 (38) 


^This assumption is made for ease of exposition. If [/ 3 ,/i] is a singular matrix, then 
further on we must use a generalized Schur decomposition for which derivations remain 
valid, but become more complicated. 
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Notice that = 0. As in the case of rational expectations models 

with constant parameters it is convenient to transform (37) using the spectral 
property of L. Namely, the matrix L is transformed into a block-diagonal 
one 


L = ZPZ-^, 


(39) 


where 


P 


A 0 
0 B ’ 


(40) 


where A and B are matrices with eigenvalues larger and smaller than one (in 
modulus), respectively; and Z is an invertible matrix®. This can be done, for 
example, by initially transforming L in a simple Schur form L = ZiLiZ^^, 
where Zi is a unitary matrix, Li is an upper triangular Schur form with 
the eigenvalues along the diagonal. We then transform the matrix Li in the 
block-diagonal Schur factorization Li = Z 2 PZ 2 ^, where Z 2 is an invertible 
matrix and P is block-diagonal and each diagonal block is a quasi upper- 
triangular Schur matrix'. Hence the matrix Z in (39) has the form Z = Z 1 Z 2 . 
We also impose the conventional Blanchard-Kan condition (Blanchard, and 
Kahn (1980)) on the dimension of the unstable subspace, i.e., dim{B) = Uy. 

After introducing the auxiliary variables 


[st,ut]'= Z ^[xt,yt^ 


(41) 


and pre-multiplying (37) by Z we have 


EtSt+\ — Ast + Qii,tSt + Q12,tUt + 

EtUt+i = But -\- Q2i,tSt + Q22,tUt + 2tEtrit+i-, 

where [Ti^t, 4 / 2 , 4 ] = and 


Qii,t Qi2,t 
Q21,t Q22,t 


Particularly, for u = 1, we have 

Ets[li = As\ ^ ^ 


= ZMtZ-^. 

+ Ql2,tu\ ^ -h \ 

+ Q22,tu[ ^ + n2,4^f \ 


(42) 

(43) 


(44) 


simple Schur triangular factorization is also possible to be employed here, but at 
the cost of more complicated derivations. The block-diagonal structure of the matrix P 
simplihes algebra 

^The function bdschur of Matlab Control System Toolbox performs this factorization. 
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where 


— + /e,*)- 

System (42)-(43) is a linear rational expectations model with time-varying 
parameters, therefore to solve the system we cannot apply the approaches 
used in the case of models with constant parameters (Blanchard, and Kahn 
(1980); Anderson and Moor (1985); Uhlig (1999); Klein (2000); Sims 
(2001), etc.). In Subsection 5.2 we develop a method for solving this type of 
models. 


ni,i 

112,t 


5 Solving the Rational Expectations Model 
with Time-Varying Parameters 

5.1 Notation 

This subsection introduces some notation that will be necessary further on. 
By I'l denote the Euclidean norm in M”. The induced norm for a real matrix 
D is dehned by 

||T)|| = sup| Ds\. 

\s\=l 

The matrix Z in (39) can be chosen in such a way that 

||A|| < a -|- 7 < 1 and ||i?~^|| < /3 + 'j < 1, (45) 

where a and (3 are the largest eigenvalues (in modulus) of the matrices A 
and respectively, and 7 is arbitrarily small. This follows from the same 
arguments as in Hartmann (1982, §IV 9), where it is done for the Jordan 
matrix decomposition. Note also that ||H||~^ < 1 for sufficiently small 7. Let 


Bt — B + Q22,t, 

At — A + Qu^f 

(46) 

By dehnition, put 

a = sup Ptll , 

b = sup \\B^ ^11 , 

(47) 




C= sup ||Ql2,ii| , 

d= sup \\Q21,t\\- 

(48) 
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In the sequel, we assume that all the matrices t = 0,1,, are invert¬ 
ible. Note that the numbers a, b, c and d depend on the initial conditions 
(xq^\zq^^). From the dehnitions of At, A, Bt, B, Qi 2 ,t and Q 2 i,t and the 
condition = (ai,0), it follows that 

lim = 0, lim d{xf‘\zf‘'^) = 0, (49) 

t^OO t^OO 

lim a{xf'\zl^'‘) = \\A\\ < 1, lim b{xf'\ zj-^'^) = < 1. 

t^OO t^OO 

This means that c and d can be arbitrary small and 

a < 1 and 6 < 1 (50) 

by choosing (xq°\ z®) close enough to the steady state. 

5.2 Solving the transformed system (42)—(43) 

Taking into account notation (46), we can rewrite (42)-(43) in the form 

EtSt+i = AtSt -|- Qi2,tUt + i,tEtTit+i, (51) 

EtUt+i = BtUt + Q2i,tSt + ^2,tEt'r]t+i- (52) 

In this subsection we construct a bounded solution to (51)-(52) for f > 0 with 
an arbitrary initial condition sq ^ 1 ^""^ and hnd under which conditions this 
solution exists. For this purpose, we hrst start with solving a hnite-horizon 
problem with a hxed terminal condition using backward recursion. Then, we 
prove the convergence of the obtained hnite-horizon solutions to a bounded 
inhnite-horizon one as the terminal time T tends to inhnity. 

Fix a horizon T > 0. At the time T using the invertibility of Bt and 
solving Equation (52) backward, we can obtain ut as a linear function of sr, 
the terminal condition EtUt+i and the “exogenous” term '^ 2 ,tEtVt+i 


Ut — —Brp^Q2l^TST — Brp^'^2,tETVT+l + Brp^ EtUt+1- 


Proceeding further with backward recursion, we shall obtain hnite-horizon 
solutions for each t = 0,1,2,...,T. For doing this we need to dehne the 
following recurrent sequence of matrices: 


KT,T-i-l — Lj,rp_- {Q21,T-i + KT,T-iAT-i) , i — 0, 1, . . . , T, 


(63) 
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where 


LT,T-i — Bt-i + KT,T-iQl2,T-ij (54) 

with the terminal condition Kt^t+i = 0. In (53) and (54) the hrst subscript T 
dehnes the time horizon, while the second subscript dehnes all times between 
0 and T + 1. Let i = 0,1,... ,T, denote the (T — i)-time solution 

obtained by backward recursion that starts at the time T. The matrices (53) 
and (54) are needed for constructing approximate solutions by backward 
recursion. 


Proposition 5.1. Suppose that the sequence of matrices (53) and (54) exists; 
then the solution to (51)-(52) has the following representation: 


UT,T-i 


'i+l 


-KT,T-iST-i + gT,i + I Lt]i 


.T 


Et-i (ut+i) , 


^fc=l 


(55) 


where i = 0,1,..., T; and 


i+l j 

9T,i = “ n -^r,r-i+fc(^2,r-i+i + KT,T-i+j'^i,T-i+j)ET-i'nT-i+j- (56) 

j=l k=l 


For the proof see Appendix A. 

The sequence of matrices (53) exists if all matrices i = 0,1,... ,T, 

are invertible. For this we need, in addition, some boundedness condition on 
the matrices Bff-KT,T-i+iQi 2 ,T-i- From (49) the matrices BfE and Qi 2 ,T-i 
are bounded, hence this condition boils down to the boundedness of matrices 

A'T,T-i+l- 

Proposition 5.2. If for a, b, c and d from (47)-(48) the inequality 


cd < 


1 

4 



/I — ab\ ^ 
) 


(57) 


holds, then 

11II ■ II A''r,T-i+i II ■ ||<5i2, T-i II < 1, i = 0,1 ,2 ,... T. (58) 

For the proof see Appendix A. 

Proposition 5.3. If the inequality (58) holds, then the matrices Lrpyp-i, 
i = 0,1,2.... ,T, are invertible. 
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Proof. From (54) and the invertibility of Bj^-i it follows that 

LT,T-i = Bx-i (/ + Bj,f.KT,T-iQl2,T-i) ■ (59) 

The matrices Lx^x-i are invertible if and only if the matrices (/ + Bff.Kx^x-iQi 2 ,T-i) 
are invertible. From the norm property and (58) we have 

\\Bx^iKx,T-i+lQl 2 ,T-i\\ < ||i?T-ill ■ l|-^r,r-i+l|| ■ ||Ql2,T-i|| < 1- 

The invertibility of (/ + Bff-Kx,T-iQi 2 ,T-i) now follows from Golub, and 
Van Loan (1996, Lemma 2.3.3). □ 

For i = T from (55) we have 


UTfi — —Kt,oSo + gT,T + 



Eq (ut+i) ■ 


(60) 


This is a hnite-horizon solution to the rational expectations model with time- 
varying coefficients (51)-(52) and with a given initial condition sq- What is 
left is to show that the solution uxfi of the form (60) converges to some limit 
as T —)■ cx). 


Proposition 5.4. If inequality (57) holds, then the limit 


lim Kx ,■ = 

T-^oo 


K, 




for j 


0 , 1 , 2 ,... 


exists in the matrix space defined in Subsection 5.1. 

For the proof see Appendix A. 

Proposition 5.5. If inequality (58) holds, then 

T+l 

n =0 (61) 

k=l 

and 

lim gx^x = 5 'oo, (62) 

T—>-oo 

where g^o is some vector in . 
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Proof. From (54) and 5.4 it follows that 


lim L 

T^OO 


T,k — 


Bk kQ12,k -^oo,/c* 


Then the limit in (61) can be represented as 

T+l T+1 

Pi n =Pi n 


k=l 


k- 


k=l 


( 63 ) 


Since iFoo,fc is bounded (it follows from formula (93) in A) and 
lim Qi 2 k = 0 , and lim BZ^ = 

k^oo ’ k^oo 

we have limk^ooLZ^i^ = B~^. Therefore, if 5 > 0 is arbitrary small, there is 
an G N such that 


\\L-^^,\\<P + 5 = p<l, (64) 

for k > N, where f3 is the largest eigenvalue (in modulus) of the matrix B~^. 
From this, the norm property and (63) we obtain 


lim 

T^oo 


T+l 

n 


k=l 


T+l 



T—>oo 

k=l 



< lim Cip'^-^ 

T-s-oo 


0 , 


where Ci is some constant. 

By (64) the products in (56) decay exponentially with the factor p as 
i —> oo. From this and the boundedness of the terms K^^k) ^ 2 ,ki 'hi ^ and 
EoPk, T G N and /c = 1, 2,..., T + 1, it follows that the series 


T+l : 

9t,t =- tt ^ 

j=l k=l 


-1 

T,k 


('h2j + KT,j'^l,j)EoPj. 


converges to some 5^00 as T —)■ oo. 


□ 


From Proposition 5.4 and Proposition 5.5 it may be concluded that as T 
tends to inhnity Equation (60) takes the form: 

Uo = —KoofiSo + Poo- (65) 

Formula (65) provides a unique bounded solution to the transformed ratio¬ 
nal expectation model with time-varying parameters (51)-(52), and may be 
treated as a policy function for this type of problems. 
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Remark 5.1. Particularly, for n = 1 we have 

T+l j 

SItIt = - E n + KTJ^lj)E«zf 

j=l k=l 
T+l j 

= - E n 

j=l k=l 

Taking into account (8), we get g^\ = 0. 

Remark 5.2. The details of derivations for the solution of time-varying ra¬ 
tional expectations model corresponding to the first order approximation of 
the system (15) and (16) are carried out in Appendix B, where we also derive 
the moving-average representation for and yl-^^. Having this representa¬ 
tion it is not hard to compute all guadratic terms in (17). 

Remark 5.3. If c = 0 or d = 0 (or both) in the ineguality (57), i.e. one of 
the variables St or Ut (or both) is exogenous to the other, then (57) is always 
valid under the conditions (50). 

Remark 5.4. The ineguality (57) is a sufficient condition for the existence 
of the solution in the form (65), and can be weakened. For the representation 
(65) we need only the invertibility of matrices defined in (54) . 


( 66 ) 


5.3 Initial conditions. 


It remains to find the initial condition for a stable solution to the system 
(51)~(52) corresponding the initial condition (16). Recall that we deal with 
the n-order problem (15)~(16), and we now put the superscript (n) back in 
notation. From (41) and (65) we have 


*0 

T^in) An) I (n) 
■-'^oo,o'’o T" yoo 


= Z 


-1 


(n) 

Vo 


where Z ^ is a matrix that is involved in the block-diagonal factorization 
(39) and has the following block-decomposition: 
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Hence 


(n) _ 712 ( 77 ) 

Vo y 

(67) 

_/y(77) in) , (n) _ '722 in) 

^ 00 , 0^0 ^ i/oo — ^ Uo ■ 

(68) 

Substituting (67) into (68) , we get 


J/S"> = (Zn + 

(69) 


The vector 0) is the initial condition corresponding to a bounded so¬ 
lution to (15) for t > 0, hence formula (69) determines the solution to the 
original rational expectations model with time-varying parameters. In other 
words, is a policy function for the rational expectations model with 
time-varying parameters at the point = 0. Particularly, for n = 1 from 

(69) and taking into account = 0 we have = 0. The condition of 

the invertibility of matrix corresponds to Proposition 1 of 

Blanchard, and Kahn (1980). 

5.4 Expected dynamics. Restoring the original vari¬ 
ables and 

To compute the expected dynamics (impulse response function) it is more 
convenient to work with auxiliary variables and then to restore the 
original variables and Substituting (65) for Ut in (51) and taking 
expectations at t = 0 gives 

BosS = (A, - + 'tijBoi)!?,. (70) 

From (67) we can compute the initial condition Sg”^ for (70). Knowing the 
initial value Sg”^ allows us to obtain the whole trajectory of the solution to 

(70) , i.e. EoStjt G K. The expected dynamics of EoUt can easily be obtained 
from (65) 

EoUt = -K^^lEoSt + Qoo- (71) 

Then the expected dynamics of the original variables is restored by 

Eox\ ^ = ZiiEosj. ^ -h Zi2EquI \ 

EoVt ^ = Z2iEos[ ^ -|- Z22Equ[ \ 
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(72) 

(73) 


where i = 1,2, j = 1,2, are blocks of the block-decomposition of the 
matrix Z. From (70) and (71) it follows that the process {st,Ut) is stable, as 
At —?• A, but A is a stable matrix, Qi2,t —t 0, and are bounded matrices 
for t > O.From this and (72) and (73) it may be concluded that the process 
is also stable. To sum up, under the assumption that the solutions 
of lower order than n are already computed in the same manner as for the 
nth order, we hnd the stable solution to the original model (1) in the form 

n 

E«y, = Y^CT'Eayf\ 

i = 0 
n 

EoXt = y^^a'EpXtK 

i=0 

6 Approximate solution: an Asset Pricing Model 

To illustrate how the presented method works we apply it to the nonlinear 
asset pricing model considered above. The simplicity of the model allows 
us to derive all approximations in the analytical form. We begin with the 
hrst order approximation determined by Equations (30) and(31) under the 
assumption that the deterministic solution and are known for t ^ 0 
and yf'^ satishes (33). Rewriting (30) for t = T and taking into account that 
-E'ra^r+i ~ gives 

y!^'^ = exp(6'a:f’|j/36' (l + exp(6'a:f’|j/3Er2/T+i, (74) 

Similarly to (74) for f = T — 1 we have 

2/r-i ~ exp{6x'^^)/36 -I- y^^'^ P^t-i + (75) 

Substituting in the last equation (74) for y!^'^ and taking into account that 
Et-ix^^ = px^^_i , we obtain 

y?-i = [0/3pexp{ex‘f^) (^1 -h yP'^ + 9{/3pY exp{e{x^p + a:®!)) (^1 + yP_^t) 

-F/?^exp(6'(a;5?^ -f- a;5?|j)ET-i|/r+i)- 
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Continuing further in the same way, for f = T — fc + lwe have 


X 


^ (i + yit) exp i 9 

.i=l \ j=l 


( 0 ) 

t+j 


xP + 


V 

= -KT,t 


( 76 ) 


(3^ exp 9 




i=i 






If the moment T tends to oo, then the following solution for is valid: 


^ exp [ 9 

i=i 


2=1 


X 


( 0 ) 

t+j 


x[^'^ = -Koc,txi^^. (77) 


=-K, 


oo,t 


Note that = 0, hence = 0. 

We now turn to the second order approximation. Equation (32) is also a 
linear forward-looking equation with time varying deterministic coefficients 
and can solved by the backward induction. Indeed, rewriting (32) for t = T 
yields 


WT 2^ yj- I yx +1 J J 

+9/3 exp{9xP_^^)Et j + /3 exp{9 x^P_^^)Et (vtIi) ■ 

Substituting (77) for y^+i in (32) and collecting the terms with Et 
yields 


y^^ = ^/3 exp{9x^^^^)9 




+/3 exp{9xP_^^)ET (^2/rii 
Substituting T — 1 for T in (79) gives 


277oo,r+i Et 


Vt-i — 1/3 exp{9x^^P)9 
+/3exp{9x'i^'^)ET-i (^yT^ 


'e (l + yf 


' ‘3‘K^ t 


Et 


(79) 


(80) 
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Inserting from (79) into (80), we have 

0(1 + 2 /®) -277oo,t 


exp(0a:®)0 


VtU 


0 [0 (1 + 2 /® 1 ) - 277oo,r+i 


, 2-/3^ exp (0 (x® + 0 :^+ 1 ) ) ^ ^ (1 + 0' 

+/32 exp (0 (x^?^ + a:® 1 )) • 

For t = T — k + 1 we have 


Et (4^^' 


Et-, (iS-li) ' 


(81) 


Vi 


'f ^ /^®xp 1^0 ^ 0(1 + 2/4) - 2iFoo,t+i ^ 

i=i \ j=i / 


(82) 


If the 


+/3^exp |^^^a;® j ^4+i- 

(2) 

moment T tends to cx), then the following solntion for y] is valid: 

«(1 + +i) - 2/1-00,,+) (83) 

provides the second term of the policy 




-i 00 

i=l 

At the time f = 0 eqnation (83) 
fnnction series expansion 


+’ = i»f;,3‘exph [»(1 + +>)-2/1 oo,.]Bo(+F (84) 

i=l \ j=l ) 

The expectation term in the last eq^uation can be obtained by nsing the 
moving-average representation for x®. Indeed, from (31) and (27) we have 

x\'^ = Ei -\- p£i-i + ... + p* 

Since the seqnence of innovations Sj, f > 0, is independent it follows that 

_i + ... + p* ^£i) 

(85) 


3 =Et{£i + p£, 


1 + p^ + ■ ■ ■ + p‘ 


,20-1) = _ _P 


1 — p2 
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The sum in exponential in (83) can be obtained from (25) 


+ ■ ■ ■ + xf'^ = X + p(xo°^ — x) + x + p^(a;® — x) 
+x + p^(xq^^ — x) = ix + ^ (a^o°^ — x). 

Finally, inserting (85) and (86) into (83) gives 



^2i 


i=l ^ 


ix + b(xQ°^ 





In computation we need to use a finite terminal time T + 1. Despite the 
fact that the method converges for any terminal condition the most 

reasonable choice of the terminal condition is the second order term in the 
expansion of the stochastic steady state in a series of powers of a. To sum¬ 
marize, we hnd the policy function approximation in the form 

y(xo) = h(xo) = I/® + 

Note that both and are functions of xq. From (83) and using (85), we 
can get the expected dynamics (in other words, impulse response function) 
Eopt. The solutions for the higher orders y[^\x), n > 2, can be obtained in 
much the same way as for yf‘\x). 


6.1 Comparison with the local perturbation 

This subsection compares the policy functions of the second order of the 
presented method with the local Taylor series expansions of orders two and 
six(Schmitt-Grohe, and Uribe , 2004). 

The parameterization follows Collard, and Juillard (2001), where the 
benchmark parameterization is chosen as in Mehra, and Prescott (1985). 
We therefore set the mean of the rate of growth of dividend to x = 0.0179, 
the volatility of the innovations to a = 0.015, the parameter 9 to —1.5 and (3 
to 0.95. For illustrative purpose, we choose the highly persistent exogenous 
process with p = 0.9 as in Collard, and Juillard (2001). 

Fig 1 illustrates the exact policy function with the approximate ones 
constructed by the semi-global method and the local Taylor series expansions. 
This figure is drawn over the interval x^ G [x — A ■ era,, x -|- A ■ cXa;], where ax 
is the unconditional volatility of the process Xt and A = 5. Fig 1 shows that 
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Figure 1: Comparison of the approximations of the policy function. The 
level of the deterministic steady state of yt is denoted by ybar 
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the semi-global approximation has the same accuracy as the sixth order of 
the Taylor series local expansion at the left endpoint of the interval under 
consideration. However, at the right endpoint of the interval the semi-global 
solution is much more accurate than the sixth order of the Taylor series 
expansion. Actually, the semi-global approximation is indistinguishable from 
the true solution in this domain. The second order of the Taylor series 
expansion is much less accurate globally than both the sixth order of the 
Taylor series expansion and semi-global solution. 

From Fig 1 one can also see another undesirable property of the the local 
Taylor series expansion, namely this method may provide impulse response 
functions with wrong signs. Indeed, the steady state value of ytisy = 12.3. 
After a big positive shock the true impulse response function is negative (the 
policy function values are below the steady state), whereas the impulse re¬ 
sponse function implied by the second order of the local perturbation method 
is positive (the approximate policy function is above the steady state). The 
sixth order approximation of the local perturbation method has the right 
sign of impulse response, but wrong shape, which is U-shaped instead of be¬ 
ing monotonically increasing. In contrast, the semi-global method, as just 
mentioned, provides almost exact impulse response. 

7 Conclusion 

This study proposes an approach based on a perturbation around a deter¬ 
ministic path for constructing global approximate solutions to DSGE models. 
Under the assumption that the deterministic solution to the model is already 
found, the approach reduces the problem to solving recursively a set of lin¬ 
ear rational expectations models with deterministic time-varying parameters 
and the same homogeneous part. The paper also proposes a method to solve 
linear rational expectations models with deterministic time-varying param¬ 
eters. The conditions under which the solutions exist are found; all results 
are obtained for DSGE models in general form and proved rigorously. 

The paper illustrates the algorithm for the second order of approximation 
using an nonlinear asset pricing model by Burnside (1998) and compares it 
with the local Taylor series expansion. The second order approximation of 
the semi-global method provide more accurate solution than sixth order of 
the Taylor series expansion around the deterministic steady state. 

The approach is applicable to Markov-switching DSGE models in the form 
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proposed by Foerster, Rubio-Ramres, Waggoner and Zha (2013), where the 
vector of Markov-switching parameters that would influence the steady state 
is scaled by a small factor. Actually, under the conditions of ’’smallness” 
of a scaling parameter and existence of higher order moments for stochastic 
terms, all derivations of Sections 3, 4 and 5 hold irrespective of probability 
distribution functions for these stochastic terms. 


Appendix A Proofs for Section 5 

PROOF OF PROPOSITION 5.1: The proof is by induction on i. Suppose 
that i = 0. For the time T from (52) we have 


Ej^ut+i — Bj^ut + Q2i,tSt + ^ 2 , rh / Thr + i - 


As Bt is invertible, we have 



where Kt^t = B^^Q2i^t] 5'r,o = — -B^^d^ 2 ,T-£'rhT-i-i and LTp\j, = BT^^. From 
(53), (54) and (56) it follows that the inductive assumption is proved for 
i = 0. Assuming that (55) holds for z > 0, we will prove it for z -|- 1. To this 
end, consider Equation (52) for the time f = T — z — 1. As the matrix Bx-i 
is invertible, we obtain 



Substituting the induction assumption (55) for UT,T-i yields 



Substituting (51) for Er-i-i(sr-i) and using the law of iterated expectations 
gives 



't\ 


2,T-iET-i-iVT-i + BrpEgx^i 
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Collecting the terms with UT;T-i-ii Sr-i-i and rjT-i-, we get 

(/ + Brp^-KT,T-iQl 2 ,T-i) UT,T-i-l = —Brp\ \{Q 21 ,T-i + KT,T-iAT-i)sT-i-l 
+ (^ 2 ,T-i + KT^T-i^l,T-i)ET-i-l'llT-i + 9 T,i + (^ 111=1 (^T+l) ] 

Snppose for the moment that the matrix 

Zt^T-i = I + Brp^-KT,T-iQl 2 ,T-i 

is invertible. Pre-multiplying the last eqnation by we obtain 

UT,T-i-l = —[(Q21,T-i + KT^T-i^T-i)sT-i-l 
+ (^2 ,r-j + Kt^T- 1^ l,T-i)ET-i-lf]T-i + gT,i 

+ (^11^=1 ^T,T-i+k^ Ex-i-l (mt+i) ] • 

Note that Lt^t-i = Bt-iZ^^t-u then using the dehnition of K^^T-i-i (53), 
we see that 


UT,T-i-l — —Kx^T-i-lST-i-l 

—Lj}rp_^ (d' 2 ,T-i + KT^T-i^l,T-i) ET-i-lTfT-i (^ 87 ) 

+Lrprp_^gT^i + Lrprp_ - ^nfc=i -^r,r-i+A:) Erp-i-i (mt+i) • 

Using the dehnition of gT,i and Lx-i^T-i+j ((54) and (56)), we deduce that 
9T,i+i = —Lrp\p_^ ( 4 / 2 ,T-i + -^r,T-i4/i,T-i) Ex-i-igr-i + Lrplp_^grp^. (88) 
From (87) and (88) it follows that 


'i+2 


UT,T-i-l — — ^Tg^-i-lST-i-l + 9T,i+l + I -f"'r,T-i-l+fc 1 Et-i-I (mt+i) • 


\k=l 


□ 


PROOF OF PROPOSITION 5.2: We begin by rewriting (53) as 

{Br-i + Kt^T-iQ 12, T-i) 77T,T-(i+l) = (Q 21 ,T-i + ET,T-i^T-i) ■ 

Rearranging terms, we have 

-^r,r-(i+i) = -Bt-j ■ (Q2i,T-i + Kt^t-iAt-i) 

— BrpEKT^T-iQ 12,T-iKT,T-{i+l) ■ 


(89) 
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Taking the norms and using the norm properties gives 


||-^r,r-(i+i)|| ^ ll-Br-ill ' IIQ2i,T-i|| + II-Bt-J ' ll-^T,r-i|| ■ ||^T-i 

+ II -By II ■ IliCr.r-jll ■ ||Qi2,T-i|| ■ ||-ft'r,r-(j+i)||- 

Rearranging terms, we get 

II jy II ^ Il-B^ijl ■ ||Q21,T-i|| + ll-Bylj ■ ||-RT,r-j|| ■ ||^T-i|| 


( 90 ) 


Inequality (90) is a difference inequality with respect to ||R''r,T-i|| > ■* = 0; 1; ■ ■ ■ > 
and with the time-varying coefficients ||y4T-i||, ||-Byij||, ||Qi 2 ,T-i|| and ||Q 2 i,T-j||- 
In (90) we assume that 

l-||RylJ-||irT,r-*||-||Qi2,T-.||^0. 

This is obviously true if ||i^'r,r-i|| = 0. We shall show that if the ini¬ 
tial condition ||iCy^y+i|| = 0, then (l - ||-By:':J| ■ ||-Rr,T-i|| ■ ||Qi 2 ,T-i||) > 0, 
i = 1, 2,..., T. Indeed, consider the difference equation: 


bd -|- basi 

(1 - bcsi)' 


(91) 


Lemma A.l. If inequality (57) holds, then the difference equation (91) has 
two fixed points 


2bd 


s^ = 


S2 = 


1 — ba + ■\/(l — bafi — Ab'^cd ’ 

1 — 6a -I- -^(1 — bafi — Ah^cd 
2bc 


(92) 


where is a stable fixed point whereas is an unstable one. Moreover, 
under the initial condition Sq = 0 the solution Si,i = 1,2,is an increasing 
sequence and converges to 


The lemma can be proved by direct calculation. From (48)-(47) the 
values a, b, c and d majorize ||Ar-j|| , ||-Byij||, \\Qi2,T-i\\ and ||Q2i,T-i||) 
respectively. If we consider Equation (88) and inequality (91) as initial value 
problems with the initial conditions ||iFyy+i|| = 0 and Sq = 0, then their 
solutions obviously satisfy the inequality ||iFyy_j|| ^ Sj+i, i = 1,2, ...,T. 
In other words, ||i7y_y_j|| is majorized by s,. From the last inequality and 
Lemma A.l it may be concluded that 


At.T-iII ^ 


1 ) 


i = 0,1,2,..., T, TeN. 


(93) 








From (92), (93) and (48) it follows that 


■ ||-ft'T,r-i|| ■ ||<5l2,T-i|| ^ 


2h‘^dc 


1 — ha + \/ 0- — haY — Ab'^cd 


(94) 


From (57) we see that 2b'^dc < (1 — a6)^/2. Substituting this inequality into 
(94) gives 




|-h"T,r-j|| ■ ||Ql2,T-i|l ^ 


(1 — bay 


2(1 -ba+ 7(1 - bay - Ah‘^cd) 


< 


(1 — ha) 1 — ha 
2(1 - ha) ~ 2 


(95) 


< 1 , 


where the last inequality follows from (50). □ 

PROOF OF PROPOSITION 5.4: The assertion of the proposition is true 
if there exist constants M and r such that 0 < r < 1 and for T G N 


WKt^j - Kt+ij\\ ^ j = 0,1,2,.... (96) 

Note now that Krpj [K-r+ij) is a solution to the matrix difference equation 
(53) at i = T — j (i = T + 1 — j) with the initial condition Kt^t+i = 0 
(R''r+i,T +2 = 0). Subtracting (89) for KT,T-(i+i) from that for 77^+1,T-p+i), 
we have 

KT,T-(i+i) — 77r+i,T-(i+i) = Brp\[KT;r-i) — 77T+i,T-i)Rr-* 

— Brp)_^KT,T-i)Ql2,T-iKT,T-{i+l) + Brp^^KT+l,T-iQl2,T-iKT+l,T-{i+l)- 

Adding and subtracting By^- ■ Kt^t-i • Qi 2 ,T-i ■ KT+i,T-{i+i) in fhe right hand 
side gives 

A'T,r-(i+i) — KT+i,T-{i+i) = Brp\{KT,T-i) — 77r+i,T-i)^T-i 
— Brp^i ■ Kt^t-i • Qi2,T-i{.KT,T-{i+i) — A'r+i,r-(i+i)) 

— Brp\{KT^T-i — KT+l,T-i)Ql2,T-i ■ A't+I.T-P+I)- 

Rearranging terms yields 

(/ + Brp\KT;T-iQl2,T-i){KT,T-{i+l) — KT+l,T-(i+l)) 

= Brj.\[KT,T-i — KT+l,T-i)A.T-i 

— Bj.\{KT^T-i — KT+l,T-i)Ql2,T-iKT+l,T-{i+l)- 
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From Proposition 5.3 it follows that the matrix 

ZT,T-i = {I + Bj,\KT^T-iQ 12,T-i) 

is invertible, then pre-multiplying the last equation by this matrix yields 

-^r,T-(i+i) — -^r+i,T-(i+i) = Zrp}j._.(Brpy{KT,T-i — KT+i,T-i)^T-i 
— Bj,\[KT,T-i) — KT+l,T-i)Ql2,T-iKT+l,T-(i+l)) ■ 

Taking the norms, using the norm property and the triangle inequality, we 
get 


||-^r,r-(i+i) — KT+i,T-{i+i)'^ 

^ ||^r,T-J ■ (ll-^T-ill ■ ll-^T,r-i — -^r+i,T-i|| • ||^r-i|| (97) 

+ II-Bt-J ■ ll-^r,r-i) — KT+i,T-i\\ ■ ||Qi 2 ,r-i|| ■ ||i^r+i,T-(j+i)||)- 

From (47) and (95) we have 


|-^T,T-(i+l) — KT+l,T-(i+l)'^ 

^ (^ab-\ - - -^ ||^T,T-ill ■ \\B^T,T-i - KT+l,T-i 

= —^—II^T,T-ill ■ \\B^T,T-i - i^r+i,r-i||- 


(98) 


From the norm property and Golub, and Van Loan (1996, Lemma 2.3.3) we 
get the estimate 


|^T,r-ill ~ 11(1 + Bj,_-KT,T-iQ 12,T-i) II ^ 


1 — ||i?2,_^iFT,T-iQl2,T-i| 




1 - ■ ||iFT,r-i|| • ||Qi2,T-i| 


By (95), we have 

liz-^ II , ^ ^ 

II T,T-^\\ < 1 + ba 

Substituting the last inequality into (98) gives 


IIB't.t-p+i) — 7LT+i,r-(i+i)|| < ||7^T,T-i — 7Lr+i,T-i||- (99) 
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Using (102) successively for i = —1, 0,1,..., T — 1, and taking into account 
Kt,t+i = 0 and Kt+i,t+i = -Bt+ 2 Q 2 i,t +2 results in 


\Ktj — Kt+ij\\ < ||-ft'r,T+i — Kt+i,t+i\\ — ||-BT+2Q2i,r+2 

^ 11-8^+211 ■ IIQ21,T+2|| ^ b\\Q21,T+2\\, j = 0 , 1 , 2 ,- 


( 100 ) 


Recall that Q 2 i,t depends on the solution to the deterministic problem (10), 

i.e. 

n - n 

V21,r — V21 I ‘^T ’ "^T+l’ 

From Hartmann (1982, Corollary 5.1) and differentiability of Q 21 with re¬ 
spect to the state variables it follows that 


||Q2i,r|| ^ C{a + 9)'^, 


( 101 ) 


where a is the largest eigenvalue modulus of the matrix A from (40), C 
is some constant and 6 is arbitrary small positive number. In fact, a + 6 
determines the speed of convergence for the deterministic solution to the 
steady state. Inserting (101) into (102), we can conclude 

\\Ktj — Kt+ij\\ < bC{a + 6*)^^^, j = 0,1, 2,... (102) 

Denoting M = bC{a -f 6) and r = a + 6 we hnally obtain (96). □ 


Appendix B The First Order System 

For n = 1 we have 


r 1 


r _di) 1 

So 

Uq 

= ZiZ-^ 

-1 

IS ^ 

_1 


From (43) for the time T we have 

= — Brp^^Q 21 ^T+lS^^ — i?2’+in2,t+l^T^ + 

Denoting Kt,t = i^^+i<52i,r+i and Rt = B:^l^Il 2 ,t+i gives 


— Rtz^^ (103) 
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For T — 1 we have 


U^T-l = -Brp^Q21,TSTli - + Bj}E t-iU^t . 


.( 1 ) 


( 1 ) 


,( 1 ) 


(104) 


Taking conditional expectations at the time T — 1 from both side (103) and 
inserting (2) we get 


Inserting (105) into (104) gives 

Mr-1 = --Br^Q2i,rSrli - B:p^Ii2,tZT-i 
+ Bj,^Et-i{—Kt,tEt-iSt — RtEz^t-^ + -Br+i-^rMr+i)- 

Inserting now Et-iSt into (106) from (42) yields 


41) 


( 1 ) 


,(i) 


“t-i ~ 


— Brp^Q2i^TST-i ~ Brp^Il2^tZT-i 


(105) 


(106) 


+ -ft"r,T(^T-iSrli + Qi2,TMrli + ni^r^rli) 

ReshnfHing terms, we have 

(/ + Bj,^ Kt^tQ 12,t)u^'P-i = —Bj.^{Q21,T + Bj,^Kt,T^T-i)St^-i 
— Brp^{jl2,t + Et,tB\,T + RtE)Z^\^ + Brp^ Bj.]^^ETU^rPj^~^. 

Multiplying (107) by (J + B^^Kt,tQi 2 ,t)~^ yields 

Mrli = ~{I + Bj,^Kt,tQi2,t) ^Brp^{Q2l,T + -BT^-^Tr^T-i)sT-i 
— (J + Bj,^K t,tQ 12,t) ^Brj,^{I\.2,t + Kj',tBi,t + RtE)z^'!_~^^ 

+ (/ + S^'iFT,TQi2,r)"'S^'5^+i^rMi^'|i. 


(107) 


or 


Mrli = —{Bt + Kt,tQi2,t) ^{Q2i,t + Bj,^K t,t^t-i) s^-i 
— [Brp + Kt,tQi2,t) ^(n2,l + Kt,tB^1,T + RtE)z^\^ 

+ (Bt + Kt,tQi2,t) ^-Br+i-^rMr+i- 
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Denoting Lt^t-i = {Bt-i + Kt,tQi2,t-i), we obtain 



Urp—X Lrprp_^{Cj2X,T Bt ,T Aj- —x) S rp 

— Lrj}rp_^{l[2^t + Kt,tAx^T + RtA)Z^-X 

Denoting 

Kt,T-1 = -^^r,T-l(Q21,T + Bj.^Kt,T^T-i) 

and 

Rt-x = Lrp]j,_-^(n.2^t + Kt^tAx^t + RtA), 

we have 



Urp _^ _ \ *b Lrprp _ -^Lrp rpE'J^Urp_^-^. 

Following the same derivation as in A for the proof of Proposition 5.1, we 
obtain the following representation: 

(108) 

where Rt can be computed by backward recursion 

Rt = (112^4 + Kj’^tB.x,t + -Rt+iA) 

Inserting (42) into (108) gives 

EtS^^i = As[ ^ + Qii^t+xs\ ^ + Qx2,t+x{—KT,ts[ ^ — Rt^t ^ 

After reshuffling we get 

F'tSt+i = (At+i — Qx 2 ,t+xET,t)s\ ^ + {—Qx 2 ,t+xRt + \ 

Denoting A* = At+x - Qx 2 ,t+xKT,t and Pi = -Qx 2 ,t+xRt + Hi^i+i, we have 

Aisli\=Aisf)+Pi4') (109) 

It is easy to see that 

D+i = ■Z’<hi+i/5,t+iD+i- 


( 1 ) 


u 


t+x 


Ai 


*t+i 

( 1 ) 


u 


t+x 
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From (109) it follows that 


{EtSt+i — St+i) + St+i — AjSt + Pi^t, 


thus, we obtain 

Si+i = AjSt + ¥tZt + 


(110) 

0 , then for f = 1 


Recall now that the initial conditions are = 0 and 
from (110) we have 




where the coefficients can be obtained by forward recursion inf = 1,2,...,T 
and backward recursion ini = 0 ,l,...,f — 1 


= At_i7t_i,t_i + P, 


t-i, 



7t,i = At_i7t_i,i + Pt_iA* ^ 


Indeed, inserting (111) into (110) and taking into account Zt = et+Eet-i+- ■ -pA* ^ei, 
we obtain 

Si+i = + 7t,t-iP-i H-1- 74,2^2 + 74,Ri) 

+Pt(£4 + Est-i + • • ■ + A* ^^i) + 

Collecting terms with Sj gives 

St+i = Pi,4^4+1 + (At7i^t + Pt)et 
+ (Ai7t,t_i + PtA)e:t_i H-h {Kat,i + PtA*"^)£i. 

Thus, for each t we compute 7 *^*, starting with the hrst index t = 1, then 
decreasing the index i = t,t — 1,... ,1 and using at each step 74 - 1 , 4 . For 
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the variable we also have a moving-average representation. Inserting the 
moving-average representation of the process and (111) in (108), we have 

Ut ^ + • • • + 7t,i^i) ~ Rti^t + R^t-1 + ■ ■ ■ + A* ^£i), (113) 

or in the shorter form 

ui ^ + • ■ ■ + ^ 4 , 2^2 + St, 1^1, ( 114 ) 

where St,i = -KT,t'yt,i - 

Taking into account that = Ziis[^^+Zi 2 u[^^ andyl^^ = Z 2 is[^^+Z 22 U^i^\ 
we get the moving-average representation for original variables 

+ ■ ■ ■ + 

Vt^'^ = P^,t^t + Pt,t-l^t-l + ■ ■ ■ + Pt,2^2 + 
where pl^ = Zii'^t,i + ^ 12 ^*,* and pl^ = Z 2 i'ft,i + ^225*,*- 
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